Time-integrated supernova neutrino flux from 

a nearby cluster 



Abstract 

The rate of gravitational collapse (type II) supernovae in our galaxy is uncertain 
by a factor of three or more. One way to determine the galactic supernova rate is 
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^ ■ through a radiogeochemical experiment (for example, the molybdenum-technetium 

. experiment) that would integrate the neutrino flux over several million years. While 

\ such a measurement is designed to integrate the flux over the entire galaxy, nearby 

00 ' star-forming regions could skew the results. We model the flux from a recently 

. identified such region, the Scorpius-Centaurus OB association, and compare with 

\ the flux from the rest of the galaxy. 
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. 1 Introduction 



Type II supernovae from the gravitational collapse of massive, short-lived stars 
are the engines that drive the chemical evolution of galaxies. While the basic 
plot of a gravitational collapse supernova (SN) is known[l], the details are 
murky: it is difficult to obtain a successful explosion in a computer model, 
and it is now generally believed that one must treat the turbulent hydrody- 
namics and/or the neutrino transport with great care in order to produce an 
outward explosion. The hot, neutrino-driven neutron wind in the aftermath of 
a gravitational collapse SN is a prime candidate for r-process nucleosynthesis, 
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but even that chapter is not finished, inasmuch as it depends sensitively on 
entropy, lepton fractions, and neutrino physics[2]. 

The famous supernova 1987A [3] in the Large Magellenic Cloud outside our 
galaxy (distance ~ 50 kpc) provided us with slightly less than two dozen neu- 
trino events; with the current generation of neutrino observatories, should a 
gravitational collapse SN occm^ within our galaxy, we will harvest a wealth of 
information to help diagnose the physics of gravitational collapse and subse- 
quent explosion (and perhaps neutrinos themselves). 

Unfortunately, at least for science, galactic SN events are rare. Various ar- 
guments generally suggest the rate of Type 11 galactic SNc to be between 1 
and 3 per century [4], although some earlier estimates were as high as 10 per 
century [5]. A galactic SN has not been visually observed for nearly four cen- 
turies, although we can only visually observe ~ 1/6 of our galaxy; in the 18-|- 
years since SN 1987A, neutrino observatories have failed to detect any galactic 
SN, which should produce hundreds or thousands of events. Knowledge of the 
galactic SN rate would help plan for neutrino observatories. 

Because the interval between galactic SNe is decades, an alternate approach 
to determining the rate is a radiogeochemical experiment [6, 7]. A radiogeo- 
chemical experiment is a generalization of radiochemical experiments such as 
the chlorine- argon and gallium-germanium experiments [8], but instead of a 
human- made target, the parent nuclei are in a naturally buried ore body and 
the daughter is a long-lived but very rare isotope. The daughter must be pro- 
duced primarily through neutrino interactions (hence one must be able to limit 
production through fission, cosmic rays, etc.), with a half-life short enough so 
no primordial atoms remain but long enough to "collect" sufficient data. As 
discussed below, and elsewhere, it is possible to identify plausible scenarios 
that will integrate the flux of neutrinos from galactic SNe. 

Gravitational collapse SNe occur only in massive (M > 8Msoiar) stars with 
short lives (typically 10-50 My), so they are associated with star-forming re- 
gions. If one uses a smooth, exponential parameterization for the disk portion 
of the Milky Way[9] one finds a flux equivalent to a point source 4.6 kpc dis- 
tant from Earth[7]. But this ignores nearby variations in the SN distribution, 
which, weighted by the inverse of the square of the distance, could have a 
disproportionate effect. 

In this paper we consider the effect of a nearby star-forming region, the 
Scorpius-Centaurus OB association (Sco-Cen), which is thought to have had 20 
SNs over a period of roughly 10 My approximately 100-130 pc from Earth[10]. 
In order to understand the effect of the SN neutrino fiux from this star-forming 
region, we consider the time-dependence of the distance to the cluster, through 
proper motion, and the finite width of the clusters. The focus of this paper 
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is on the relative contributions from nearby clusters, Sco-Cen in particular, 
relative to the neutrino flux from SNe throughout the galaxy. 



2 The Mo-Tc experiment 

The example we consider is the molybdenum-technetium (Mo-Tc) experiment [6 , 1 1] ; 
other potential radiogeochemical experiments are discussed in Ref. [7]. The 
Mo-Tc experiment was originally proposed to measure the long-term solar 
neutrino flux, in order to test the hypothesis that the measured low flux of 
solar neutrinos might be caused by nonstandard solar models. Such nonstan- 
dard solar models have been ruled out through helioseismology and detection 
of neutrino flavor oscillations [12]. Inverse beta decay induced by neutrino flux 
on molybdenum ore produces technetium 

^^Mo + i/^^^Tc + e-. (1) 

^^Tc has a half-life of 4.2 My, short enough so that it cannot be primordial in 
origin, but long enough to integrate over the mixing time of the solar core. 

Originally the main source of neutrinos considered were neutrinos from the 
Sun. It was later realized[7] that ^^Tc, with a half-life of 2.6 My, could serve 
as a tracer of supernova neutrinos through the reaction 

^^Mo + u Tc + n + e-. (2) 

Solar neutrinos are generally too low in energy, however, to knock-out a neu- 
tron and produce '^^Tc through (2), which has a threshold of 7.28 MeV, but 
more energetic SN neutrinos can easily do so. Table 1 gives the relevant cross- 
sections as computed in [7]. (In this paper we do not consider neutrino flavor 
oscillations; such oscillations will certainly change cross-sections[13], but focus 
only on relative contributions to supernova fluxes.) 

The larger cross- sect ions of SN neutrinos allow them to compete with solar 
neutrinos as a source of ^^Tc despite having a lower flux. The standard solar 
model flux of ^B solar neutrinos is ~ 6 x 10^ cm~^ s~^ [8] (this does not include 
neutrino flavor oscillations, which will reduced this considerably), while the 
flux of galactic SN neutrinos is much smaller, about 

X iVss X 1.3 X 10=*cm-V\ (3) 

where R is the rate of galactic gravitational collapse SNe per century, some- 
where around 1 to 3, and N^s is the number of neutrinos emitted per collapse 
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in units of 10^^, about 0.45 for electron neutrinos. (A similar study [14] looks 
at the neutrino flux from thermonuclear burning throughout the galaxy and 
from extragalactic sources.) 

This beautiful dichotomy is spoiled somewhat by the existence of ^^Mo (9.55% 
compared to 24.13% for ^^Mo), so that 

^'Mo + u,oiar^^'Tc + e- (4) 

can also contribute. The cross-section is not known because there is no experi- 
mental data on the Gamow-Teller strengths, but a reasonable scaling argument 
[7] can be used to derive the ^^Mo cross-sections in Table 1; one would need 
to subtract off the solar contribution in order to measure the galactic SN rate. 
(An attempt to mount the Mo-Tc experiment was begun by K. Wolfsberg and 
collaborators in the 1980 's, but the original experiment was never completed 
due to funding hmitations [15]. We note, however, that the Henderson molyb- 
denum mine in Colorado, the site of the original Mo-Tc experiment, is still 
active and is at the time of this writing a candidate for siting the National 
Underground Science Laboratory, and it is not impossible that the Mo-Tc 
experiment could be resurrected.) 



3 Time-dependence and nearby events 

Sco-Cen can be divided into through subgroups: Lower Centarus Crux (LCC), 
Upper Centaurus Lupus (UCL), and Upper Scorpius (US), each of whose 
distance from the Earth varied between 100 and 240 pc over the past 10 
My[16,10] due to proper motion. Each of the subgroups have a finite window 
over which SNe could have occurred: 3 My for US, 7 My for LCC, and 10 My 
for UCL. The SN rate in each subgroup is assumed to be approximately 1 per 
My[16,10]. 

In order to account for the time-dependence of the neutrino flux, we note we 
must solve 

^^-XN{t)+am, (5) 

where N{t) is the number of daughter nuclei, here ^^Tc, A is the decay constant, 
a is the reaction cross-section, and ^{t) is the time-dependent flux. From the 
data in Ref. [10] we approximated $ as a quadratic in t; this can be solved 
analytically. We assumed the rate of SNe is constant in time. 

The most useful way to present the results is to compute an effective distance 
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Tes; the effect of the timc-intcgrated flux for the cluster in motion is the same 
as a static cluster fixed at the distance Teff. 



We found that the effective distance Ves can be approximated by the following 
simple prescription: 



where T is the time before the present that SNc began in the subgroup; this 
formulation, which is simpler than solving Eq. (5), agrees with such solutions 
to within a fraction of a percent. Table II gives the average fluxes and Teff for 
each subgroup studied. The total average flux from Sco-Cen is then the sum 
of the fluxes in the middle column, 450 cm~^ s~^. 



4 Finite-size corrections and extra-near events 



We also considered the fact that the star-forming subgroups have a finite 
width, on the order of 20-30 pc. By considering a sphere of uniform density 
with a center at a distance r from Earth and with rms radius a, we find that 
the flux from the cluster is 



Monte Carlo simulations with a Gaussian distribution but the same rms radius 
agreed with this correction term. Thus the correction due to flnite width is 
relatively small. 

We further took the opportunity of our Monte Carlo simulation to investigate 
the question of extra-near SN events. Anomalous deposits of ^°Fe in deep-sea 
crust have been linked to a nearby (a few 10 pc) SN within the last few My 
[17,18]. Benitez, Maiz-Apellaniz, and Canelles [10] suggest this anomaly could 
be explained by a SN in Sco-Cen; taking a cluster at 100 pc with a Gaussian 
with an rms width of 25-30 pc, they argue for a SN as close as 40 pc from Earth 
(2a lower limit). Unfortunately this argument is based upon a one- dimensional 
Gaussian distribution and is wrong for a three-dimensional distribution. In our 
Monte Carlo simulations, we found relatively few simulation runs that placed 
a SN at about 40 pc; in fact only a fraction of a percent of any given simulation 
with 20 SNe gives an event closer than 50 pc. Increasing the rms width of the 
cluster to 40 pc still only yields a near-Earth ( < 50 pc) event only 12% of 
the time. Of course, such simulations depend sensitively on the tail of the 
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distribution, which might be non-Gaussian. Nonetheless, while Sco-Ccn is not 
ruled out as a candidate source for the ^°Fe anomaly, the statistical linkage is 
weaker than previously argued. Conversely, one could possibly use ^°Fe and 
other radioisotopes to constrain nearby events [19] and sharpen our modeling 
of the flux from such events. 



5 Conclusions 

We have considered the potential effect a nearby star-forming region, Sco-Cen, 
would have on radiogeochemical measurements of the galactic supernova flux, 
in particular the Mo-Tc experiment. Although we performed detailed calcu- 
lations, we found that both the time-dependence of the flux due to proper 
motion and the finite width of a region can be calculated through simple, 
effective prescriptions which could be applied to any cluster and any radio- 
geochemical experiment of interest. The flux from Sco-Cen is roughly the same 
as that from throughout the galaxy if the galactic SN rate is 1 per century; 
this "background" could limit the ability of a radiogeochemical experiment to 
measure the overall rate. At higher galactic rates the Sco-Cen contribution 
will be relatively small. 
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Table 1 

Neutrino capture cross-sections for various sources; taken from [7] 

Reaction v source a ( 10~^^ cm^) 

98Mo(z/,e-)98Tc Solar 8b 3 

98Mo(z/,e-)9^Tc SN 50 

98Mo(z/,e-)9^Tc+n Solar 0.3 

''^Mo(z/,e-)9'^Tc+n SN 150 

9^Mo(i/,e-)9'^Tc Solar 1.2 

97Mo(z/,e-)97Tc SN 18 



Table 2 

Subgroups in Sco-Cen 

Subgroup T (My) ($) (cm~^ s"^) r^^ (pc) 

US 3 111 154 

LCC 7 218 110 

UCL 10 120 148 
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